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Abstract. We use the method of statistical parallax to constrain the distance- 
scale zero points and analyze the kinematics of extensive samples of Galactic 
classical Cepheids, RR Lyrae type variables, and blue horizontal branch stars, 
which serve as standard candles/kinematic tracers of various Galactic popula- 
tions. We obtain three consistent estimates for the local circular velocity based 
on the mean velocities of halo RR Lyrae variables, BHB stars, and Galactic rota- 
tion curve inferred from Cepheid data with an average value of 210 ± 6 km s~^, 
which is close to the average circular velocity in the 5-40 kpc interval of Galac- 
tocentric distances inferred from BHB star data (195 ± 5 km s~^), thereby pro- 
viding further supporting evidence for the practically flat shape of the Galactic 
rotation curve beyond ~ 5 kpc from the center. The inferred distance-scale cor- 
rections imply a solar Galactocentric distance of 7.7 ± 0.4 kpc, an LMC distance 
modulus of 18.42 ± 0.06, and a Hubble constant of 73-85 km 8"^ Mpc'^. 



1. Introduction 

Establishing the kinematics of the major Galactic subsystems — the halo and 
the thin and thick disk — is a task of fundamental importance, bringing us closer 
to understanding the origin, evolution, and mass distribution in the Galaxy. In 
this work, we choose Cepheids as our kinematic tracers for the young (thin) 
Galactic disk and RR Lyrae variables as the kinematic tracers for the thick disk 
and halo. We further expand our halo kinematic-tracer base by including an 
extensive sample of blue horizontal-branch stars (BHB) spanning heliocentric 
(and Galactocentric) distances from ~ 5 to several tens of kpc. 

Here we use the method of statistic a l par allaxes in its rigorous maximum- 
likelihood version sugg ested by iMurravl (|l983l l (pp. 297-302) to estimate the 



kinematical parameters of the sample of Galactic Cepheids, RR Lyrae type vari- 
ables, and BHB stars and to calibrate the underlying photometric distance scales 
(based on the period-K-band luminosity relations for Cepheid s and RR Lyraes 
and on the BHB star distance scale proposed by ISirko et al. I ( 20041 ) . which is 
based on SDSS colors and the appropriate model atmospheres). 



2. Cepheids 

Our sample consists of 215 Cep heids with bona-fid e proper motions from the 
new reduction of Hipparcos data ( van Leeuwen |[2007l ) . accurate 7- velocities from 
Gorynya et al. (1992, 1996, 1998) and photometry from Berdnikov (1995, 1998) 
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and Berdnikov et al. (2 000). The distanc e s are based on the Cepheid ivT-band 
PL relation derived by iBerdnikov et al. I (|l996 l) and the procedure proposed 



in the latter paper. The results are listed in Table 1. Here Uq, Vq, and Wo 
are the components of the local mean solar velocity toward the Galactic cen- 
ter, in the direction of Galactic rotation, and toward the North Galactic Pole, 
respectively, relative to the sample considered, ctVr, crV^p, and aVz are the diag- 
onal components of the velocity dispersion tensor in the directions of projected 
Galactocentric radius. Galactic rotation, and toward the North Galactic Pole, 
respectively. The parameters ^Iq, A, and Q" are equal to the local angular cir- 
cular rotation velocity, Oort's constant A, and the second derivative of the local 
angular circular rotation velocity with respect to projected Galactocentric dis- 
tance, respectively, and k is the distance-scale correction factor (to be multiplied 
by the initial distances to obtain the corrected distances). 

Table 1. Kinematical and distance-scale solution for Cepheids (N=215 
stars). Here Uq, Vq, Wq, ctVr, ctV^, and aVz are in km s~^, flo and A 
are in km s^^ kpc^^, and il" is in km s~^ kpc^'^. 



Uo Vo Wo (tVr aV^ aVz VLq A U" k 

"^8:5 ^2:2 n 131) 9:3 8^9 27:5 TKT 080 Lis 
± 1.0 ± 0.9 ± 0.9 ± 0.8 ± 0.6 ± 0.9 ± 0.8 ± 0.7 ± 0.13 ± 0.06 



3. RR Lyrae type variables 



Our sample consists of 364 RR Lyrae type variables with bona-fide proper mo- 
tions , radial velocities , and phase-corrected 2MASS K-ba nd magnitud e s. Th e 
dist ances ar e based on the i^T-band PL relation derived by I Jones et al. I (|l992l '). 
See Dambid (j2009. ) for a complete description of the data set, which is available 
from the CDS. We ap plied to it a bimodal version of the statistical-parallax 
method ( Dambisll2009l ) and obtained the results that we list in Table 2. 



4. BHB stars 



Our third sample consists o f 1955 rigorously selected blue horizontal-branch 



halo stars from SDSS DR6 (lAdelman-McCarthv et al. 



ric distances based on the calibration of Sirko et al. 



.20081 ) with photomet- 
200411. radial velocities 



determined within the fr amework of SEGUE project ( Yannv et al. 1120091 ). We 



adopted this sample from lXue et al. I (12008.') . and supplemente d it with absolute 



proper motions from the SDSS DR7 ( Abazaiian et al. 20091 ) database. Here 



we use only stars sample with Rg < 40 kpc and removed a number of evi- 
dent outliers. We applied to the resulting sample a single-mode version of the 
statistical-parallax method parametrizing the variation of the velocity-ellipsoid 
parameters as a funct i on of Galactocentric distance in terms of the model of 
Sommer-Larsen et al. ( 19971 ) and obtained the results that we list in Table 3. 
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Table 2. Kinematical solution and distance-scale correction for Galactic 
field RR Lyrae variables based on the bimodal solution (364 stars with rhei < 
6.5 kpc). Here Uq, Vq, and Wq, have the same meaning as in Table 1, and 
<^Vii, (tV^, and aVg are the diagonal components of the velocity dispersion 
tensor aligned with the spherical coordinate system centered on the Galactic 
center. All these quantities are in km s^^ 



Population 


Fraction 


Uo 


Vo 


Wo 


ctVr 






k 


Halo 


0.75 


-12 


-217 


-6 


167 


86 


78 






0.03 


± 10 


± 9 


± 6 


± 9 


± 6 


± 5 


0.97 


Thick disk 


0.25 


-15 


-45 


-25 


55 


44 


30 


± 0.04 




0.03 


± 7 


± 7 


± 5 


± 6 


± 6 


± 4 





Table 3. Kinematical solution and distance-scale correction for the sample 
of BHB stars (1955 stars with r^ei < 40 kpc). Here Uq, Vq, and Wq have the 
same meaning as in Table 1 a.nd an, a~^, I, rp, and Vc are the parameters of 
the kinematical halo model of ISommer-Larsen et al. I (|1997[ ). The parameters 
Uq, Vb, Wo, ctq: cr'*', and Vc are in km s~^, and I and rp are in kpc. Note that 
Vc is the circular rotation velocity of the (assumed) flat rotation curve and 
ro is the Galactocentric distance at which the shape of the velocity ellipsoid 
changes from radial to tangential anisotropy. 





Vo 


Wo 


(Jo (J+ 


I 


ro 




k 


-14 


-225 


-2 


96 89 


3.7 


19.2 


195 


1.06 


± 4 


± 5 


± 3 


±12 ±9 


± 0.4 


± 1.6 


± 5 


± 0.03 



5. Distance scale 

In Table 4 we list the distance-scale correction factors inferred via statistical 
parallax (k) and and the average distance-scale correction factors based on sta- 
tistical and Hipparcos trigonometric parallaxes (^2) for Cepheids and RR Lyrae 
variables, as well as the corresponding LMC distance-modulus estimates. We 
also give the estimate of the solar Galactocentric distance (Ro) based on the 
IR photometry of Galactic bulge RR Lyraes. The weighted average of the 
LMC distance modulus {DMlmc=^8A2±0.0Q) implies a Hubble constant of 
73-85 km s"^ Mpc'^ 



Table 4. Distance-scale results. 



Distance indicator 


k 


k2 


DMlmc 


Ro, kpc 


Cepheids 


1.18 ± 0.06 


1.12 ± 0.03 


18.50 ± 0.07 




RR Lyraes 


0.97 ± 0.04 


0.98 ± 0.04 


18.30 ± 0.09 


7.7 ± 0.4 
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6. Circular velocity 

Our kinematical analysis yields three independent estimates for the local velocity 
of circular Galactic rotation, Vc{local): Vc{local) = \Vo{halo) — Vo{Cepheids)\, 
where Vo{halo) can be represented either by halo RR Lyraes or BHB stars (the 
corresponding two Vc{local) values are equal to 205 ± 9 and 205 ± 13 km s~^, 
respectively), or as Vc{local) = \Ro ■ QqI = 212 it 13 km s^^. Tlicy agree rather 
well with each other and their weighted average Vc{local) = 210 ± 6 km s~^ 
differs little from the global value of Vc{global) = 195 it 5 km inferred from 
an analysis of the kinematics of BHB stars out to a Galactoccntric distance of 
Rg ~ 40 kpc, thereby further reinforcing the case of the flat Galactic rotation 
curve beyond -Rg ~ 5 kpc. 
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